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I. INTRODUCTION 


In quantum gravity, we are interested in the computation of transition amplitudes between spatial geometries. 
These transition amplitudes are supposed to arise as a sum over spacetime ‘geometries’. In the classical realm what 
is meant by geometries is clear, one usually means a smooth Lorentzian metric free of singularity on a manifold of 
fixed topology, possibly satisfying some asymptotic conditions. In the quantum realm we obviously have much more 
freedom in choosing which spacetime geometries should be included in the sum. Do we authorize the geometries to 
posses some singularities for instance? [1] Do we authorize topology change? One of the outstanding question in this 
context is wether one should or not include a sum over topology in the sum over geometries. The next outstanding 
question is wether this is possible in a given framework. In this paper we are going to analyze this problem in the 
context of discrete approach to quantum gravity, i-e dynamical triangulation and spin foam models. In dynamical 
triangulation one first choose an equilateral triangulation A of space-time with regular simplices and we assign to 
each triangulation a statistical weight given by the exponential of the Regge action. In spin foam models one choose 
a triangulation of space time and assign to it certain representation labels of the Lorentz group. This representation 
labels carry information on the geometry of the simplex. The transition amplitudes are constructed as a state sum, 
using sum over representations assigned to the geometrical building blocks of the triangulation. Historically, the first 
spin foam model was proposed by Ponzano and Regge [2] for 3D Euclidian quantum gravity, and their work has been 
recently generalized to higher dimension to give birth to the spin foam approach for quantum gravity (see [3, 4] for 
an overview). 

3D gravity is a topological theory. As a consequence, the transition amplitudes built in the Ponzano-Regge model 
do not depend on the choice of the triangulation. However this is no longer the case in 4D, and one has somehow to 
get rid of the choice of the triangulation in order to restore the infinite number of degrees of freedom of 4D gravity. 
A proposal to do this is to sum the discretized amplitude over all possible triangulations of the manifold. A more 
drastic proposal is to sum over all possible triangulations regardless of the topology. 

It is however commonly believed that there is no way to give sense to the sum over triangulations including all 
topologies. The problem appears clearly in the context of dynamical triangulations [5]. In this context the amplitude 
of a given 3-dimensional triangulation A is given by 


A[A] 


1 

Sym[A] 




( 1 ) 


where Sym[A] is a symmetry factor, and (resp. Ni) denote the number of tetrahedra (resp. edges) of the 
triangulation. The parameters A, ^ can be interpreted, if one uses the Regge action, in terms of the bare cosmological 
constant A and the bare Newton constant 1/G 


exp(^). 

(2) 

A ~ exp(-A- ^), 

(3) 

where c is a geometrical constant. The sum over all triangulations is then 



(4) 

N3,Ni 



with Af {Ns, Ni) being the number of triangulations with a fixed number of tetrahedra and edges. Since /i > 1 the 
partition function is bounded from below by X^^J\f{Ns) where is the total number of triangulations with 

a fixed number of tetrahedra. This number grows factorially with the number of tetrahedra if there is no restriction 
on the topology of the triangulation. Therefore a sum over triangulations does not seem to be possible in general, 
since this factorial growth can not be killed by any choice of A, /i. The solution proposed in [5] is to restrict to a fixed 
topology (as the 3-sphere), in order to get back to a case where the number of triangulations grows exponentially with 
the number of 3-simplices. So the sum over simplicial manifold can be defined in the case of a restricted topology. 
However, since there is no general classification of the topologies in D > 2, we cannot expect to use this result to 
make sense of a sum over topologies as in 2 dimensions. One conclusion that one can draw is that one should not 
try to sum over topology when trying to solve quantum gravity. This is not completely satisfactory, since first, there 
are argument in favor of classical topology change [1], but also one would expect in a theory of quantum gravity to 
resolve this issue dynamically and not a priori. 

There is hopefully another point of view that can be taken on this problem. It is well known since Dyson [6] 
that physical quantities in interacting quantum theory are non-analytic functions of the coupling constant. This 
translates into the fact that the perturbative expansion of any physical quantities in terms of the coupling constant 



3 


is a divergent series. This opens the possibility of giving a meaning to the sum over triangulations if we can interpret 
it as a perturbative expansion of a non-perturbative quantity. 

These ideas have already been used in the context of 2D gravity, leading to the construction of the matrix models, 
which realize a sum over triangulations of 2D surfaces [7]. By considering three-indices tensors instead of matrices, 
Ambjorn et al. [5] tried to apply these ideas to obtain a dynamical triangulation model for three-dimensional quantum 
gravity. Later, Boulatov and Ooguri [8, 9] have generalized these models in the form of non local field theories over 
group manifold. The bottom line is that a space-time triangulation can naturally be interpreted as a Feynman graph 
of a Group Field Theory (GFT) [10]. Moreover, the specific gravity amplitude assigned to a given triangulation is 
exactly given by the Feynman graph evaluation of a group field theory GFT. This was shown in the context of 3D 
dynamical triangulation [5], 3D Ponzano-Regge model [8] and even 4D spin foam model [11]. It is now understood 
that this is a general feature of spin foam models [12]. 

All these results lead to the conclusion that the sum over triangulations of a state sum amplitude is in fact a sum 
over Feynman graphs of a GFT, 

“Estate sum(^) = “^GFT(r)- (5) 

A r 

Therefore it is not a surprise that the sum over triangulations is divergent, since it is a perturbative expansion, it 
should be understood as an asymptotic series. If we take for serious the fact that the sum over triangulations comes 
from a group field theory, this opens the possibility of giving a non perturbative meaning to the sum. This idea has 
been already advocated in [12] and to some extent in [5]. 

However it is well known that in general, the asymptotic series does not determine uniquely the resummation. One 
need additional non-perturbative inputs. In interacting quantum field theory, the large order behavior is determined 
by non-perturbative information contained in the action and in some case [13] one is able to reconstruct a non- 
perturbative expression from the divergent series, in a uniquely defined manner. In such cases, the divergent series 
is said to be uniquely Borel-summable. Our main result in this paper is to prove that such a procedure can be 
realized for a particular type of GFT models associated with three dimensional gravity. If the Feynman expansion 
of a GFT model is Borel-summable, this means it could be related to a non perturbative expression in a uniquely 
defined manner, giving a well-defined meaning to the equality 

( 6 ) 

r 

This result proves that a well-defined notion of a sum over all triangulations of all topologies is not unreasonable, as 
it could be achieved in particular cases. 

In the context of dynamical triangulation the coupling constant A of the group field theory has a physical meaning, it 
can be interpreted as — exp(—A), where A is the cosmological constant. Since it is negative, the potential is unbounded 
from below and it was expected that the corresponding GFT is not uniquely Borel summable in the physical regime. 
Therefore one can think that the theories allowing non perturbative resummation will not be physical (at least in 
the context of dynamical triangulations). In two dimension one can view this problem as the one preventing a non 
perturbative definition of string theory via matrix models. 

It was therefore for us a big surprise and a relief to find out that the modification of the Ambjorn et al. model [5] 
we propose is uniquely Borel summable in the physical region where A > 0 (see sect. IV). In the context of three 
dimensional Ponzano-Regge models the value of the group field theory coupling constant model is not restricted by 
physical requirement and we expect the Borel summable GFT to have physical meaning without restriction. 

In section II of this paper we first review the definitions of several group field theory models related to 3 dimensional 
gravity like the Boulatov or Ambjorn model. We present their Feynman expansions and show their relations with 
3 dimensional Ponzano-Regge models or dynamical triangulations. We also discuss their convergence properties. In 
section III we review known facts and theorems concerning non perturbative resummation of a perturbative series. 
In section IV we propose a modification of the usual Ambjorn or Boulatov models in order to deal with positive 
potential. We show that the perturbative expansion of this model can be related to the dynamical triangulation sum 
with physical value of the parameters. In the last section V we finally prove the unique Borel summability of the 
model we have defined. This give us an example of an unambiguous way to sum over all triangulations hence all 
topologies. 


II. 3D GROUP FIELD THEORY MODELS 

In this part, we review various constructions of 3D group field theory (GFT) models : the Boulatov model, leading 
to the Ponzano-Regge amplitude [8], a model on the non-commutative sphere, reproducing the 3D matrix model 
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of Ambjorn et al.[5], and some others models with improved convergence properties, analogous to those recently 
introduced for 4D spin foam models [14, 15]. 


A. The Boulatov model 

The Boulatov model is a field theory defined on the group manifold S'0(3)^, with some invariances imposed on the 
fields. The fields are asked to be real ($ = $), symmetric 

^{9a(l),9a{2),9a{3)) = (7) 

(where a denotes a permutation of three elements and |ct| its signature), and (right-) SO{3) invariant 

'^{919,929,939) = ’^{91,92,93)- ( 8 ) 

This last requirement allows us to see the Boulatov model as a field theory on the coset space SO(3)^/SO(3). The 
action is non-local and is given by 

- 5 '[^'] = ^ / dgidg2dg3 <^{gi,g2,93f 
^ Jg3 

+ 77 / dgidg2dg3dg4dg5dge ^{gi,g2,g3)^{93,95, 94 )^( 94 , 92 ,96)^(96,95,9i) (9) 

4! Job 

To compute the partition function of this model in perturbative expansion, one has to express the action in terms 
of the Fourier coefficients of the field 4). In order to define a Fourier transform on the coset space, one first considers 
the transform on SO(3)^, then requires the SO(3)-invariance. The representations of SO(3) are labelled by integer 
spin j, and a function on SO (3) can be expressed as a sum over the matrix elements of g in the spin j representation 
DLn(9) (see appendix A) 


= E dj^rOi^^ig), (10) 

j,m,n 

where dj = 2j + 1 and — j < m,n < j. The Fourier coefficient is 

^ f <i>(g)Di,^(g*), ( 11 ) 

JG 

where g* denotes complex conjugation. The Fourier expansion on SO(3)^ is then 

H9i,92,g3) = E (12) 

j,rh,,n 

where J = (ji, j 2 , J3) and Dl^^ig) = (5i)71^2.«2 (52)71^3 ,„3 (93)- The SO(3)-invariance is implemented by group 

averaging over the right diagonal action of SO(3) 

*^( 91 , 92 , 93 ) = J dg <^{gig,g2g,939)- (13) 

Expanding the RHS using decomposition (12), then performing the integration on g using the relation (A6), we obtain 
the Fourier decomposition of the field over SO(3)^/SO(3) 

H9i,92,g3) = E (14) 


where we introduced the Fourier coefficient 

n 


(15) 
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FIG. 1: Graphical representation of the propagator and potential. The black square stands for the sum over different permu¬ 
tations of the three lines, see (20). 


and Ci denotes the Wigner 3j-symbol (see appendix A). The symmetries asked on the fields are translated on the 
coefficients : 


7 mim 2 m 3 

3 (T (1) 3 (T (2) 3 (T (3) 

We can now compute the action (9) in terms of the Fourier coefficients. We insert the Fourier decomposition (14) in 
the action, and perform the integrations using (A5). The four 3j-symbols appearing in the potential term recombine 
into a 6 j-symbol. We finally obtain the action in Fourier modes 

C _ _j_ \ ' I 3mim2"»3 |2 \ ' /I 3 -mi-m 2 -m 3 <033-0357714 <-034032703 A-rngm^mi 

2 2 -^ I iiAjs I 41 <^ 1 1 jihJa jsJaji hhie, AAii 

j ji T •-je ,Tni. .rriQ 

( 18 ) 

We want to compute the partition function and correlation functions of the theory. They are formally defined in 
terms of the path integral 


f Ji j 2 js 1 
1 J 4 js J6 j 


I" _ 2 i (it-l-TOi ) ^-701 - 032 - 703 

(-i)|"Ia:; 


jmim2m3 

'^jlj2j3 


(16) 

(17) 




(19) 


As usual we compute this partition function as a perturbative expansion. The Feynman rules of such a group field 
model are very similar to those of a usual QFT. One first construct the propagator and then include the interaction 
perturbatively. The propagator is given by 


< AfA^ >= 1 y , 

j j 31 




( 20 ) 


where the sum is over permutation of three elements. The full contribution can be recast into a sum over Feynman 
graphs, whose form is dictated by the action. For each graph we compute its Feynman amplitude by labelling the lines 
of the graph by representations indices; associate a weight to the lines and the vertices; then sum over representation 
indices. In our model the graphical form of the propagator and the vertex are given in figure 1. The contributions 
of the theory at order n are obtained by constructing graphs with n vertices. We then label each closed line of the 
graph by representation indices (j, m ), associate to each vertex the weight given by the 6 j-symbol, then sum over 
representation indices. Considering a graph with n vertices and a symmetry factor sym[r], we obtain the amplitude 


A[r] 


(~^)" TT / \ 

VI i 


( 21 ) 


As the weight does not depend on the indices m on the lines, the sum over each of the m indices is trivial and equal 
for each line to the number of terms dimj) = 2ji + 1. Thus one get the amplitude 


A[r] 


(-A)" 

sym[r] 


En(27<+i)n 

fil I 


Jl J 2 J 3 1 

J4 3b je j 


( 22 ) 


By associating a triangle to each propagator and gluing them together at vertices, one generates tetrahedra (as seen in 
figure 2) and in general, a Feynman graph of the theory is interpreted as a 3D triangulation. The state-sum amplitude 
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FIG. 2: Triangulation generated by the Feynman diagrams 


associated to a triangulation A is thus 

bl e t ^ ^ 

where e and t label the edges and tetrahedron of the triangulation. This amplitude is exactly the Ponzano-Regge 
partition function associated to the triangulation A. Moreover, all the triangulations can be obtained as Feynman 
graphs of the theory, thus the sum over Feynman graphs is indeed the sum over triangulations of the Ponzano-Regge 
amplitude. 

As in standard quantum field theory, the GFT are in general faced with two kind of divergences. The first ones 
are the analog of the UV/IR divergences, and arise when computing the value of a given graph by summing over 
representation indices (momenta). The second ones are the large order divergences, arising when one try to sum all 
the Feynman graphs, in order to get a non-perturbative expression of a physical quantity. In this paper we will mainly 
concentrate on the large order divergences, and we need somehow to get rid of the UV/IR divergences. For instance, 
it is well known that the weight (22) is divergent for some triangulations, the divergences being associated with 
the vertices of the triangulation. This is a generalization of the usual fact that in standard field theory divergences 
are associated with the loop of the Feynman graphs (for more details on loop divergences in group field models see 
[14, 15].) As we will see at the end of this section, one way to avoid these divergences is to define spin foam models 
with enhanced convergence properties, in analogy with the models introduced in [14, 15]. However the Ponzano-Regge 
weight generated by the Boulatov model is not one of these enhanced models, and we feel that a full understanding of 
the notion of renormalisability for these theories should be developed. In this paper we will not dwell further on this 
issue and we will ignore the UV problems by working with a fix regulator and concentrate on the problem of large 
order divergences. 

The natural regularization to choose in this context, which respects the symmetry, is to project out modes of the 
fields carrying a spin j bigger than say N. So instead of the field $ we use the projected field 

$iv(ffi,ff2,ff3) = E (24) 

j^N fn,n 

In this case the action expressed in Fourier mode is identical to (18), except that the sum over the spins are restricted 
to be over j < N. One can naturally interpret this projection as arising from a group field theory defined on the 
non commutative sphere [16] (see [17] and the appendix B for non commutative 2-sphere S%). This leads to an 
interpretation of the regularized group field theory model in terms of a generalized matrix models, since functions on 
the non commutative sphere are matrices. Another natural way to get a regularization is to define this model with 
the q-deformation of SU(2) for q a root of the unity q = where the cut off N in the representations is related to 

the cosmological constant A by 2N -I- 1 = ttA”^/^. 


B. Group field model on the sphere and 3D matrix models 

Keeping the action (9), one can define various models by requiring different kind of invariances on the fields. For 
instance by considering a real symmetric field on SO(3)^ with a SO(2)^-invariance, 

4‘(ffi^i, 52 / 12 , 53 / 13 ) = 4>(5i, 52,53) V(/ii,/i2,/i3) G SO{2f (25) 

one obtains a group field model on the coset space [SO(3)/SO(2)]^. To implement the SO(2)^-invariance, one fixes 
an SO(2) subgroup of SO(3). In particular, we may consider SO(3) in its fundamental representation and choose a 
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unit vector v^; we then let H be the SO(2) subgroup which leaves invariant. The space of the equivalence classes, 
SO(3)/SO(2), is diffeomorphic to the 2-sphere 8“^. Thus a GFT over the group manifold [SO(3)/SO(2)]^ is viewed as 
defined over three copies of the 2-sphere S^. As in the case of the Boulatov model, one has to determine the Fourier 
transform by requiring the invariance on the Fourier expansion. If the representation is of integer spin j there is a 
unique vector invariant under the action of the SO(2)-subgroup, we denote it The SO(2)-group averaging of the 
matrices of representation gives the projector over this invariant vector 

[ dh Dl^ih) = vlvl (26) 

JS0{2) 


So if we define the Fourier coefficient 


A?" 

3 


vT 

n 


(27) 


the field is expressed as 


$(<71,52,53)= 

j,m,n 


(28) 


The reality conditions are translated on the coefficients 


3J1J2J3 _ ll lEi 3J102 J3 

V ■*-/ mi—m 2 —m 3 


The action is found to be 


(29) 


S = 


- E 

jij2j3,mi,m2.m3 


Amim2m3 

^313233 


X 

4! 


E <->F‘ 


"i A— — "^2—"^3 Am3 — m5m4 *—m4m2mQ A—mem^mi 
^3i32j3 ^j33534 '^34:j2j& ^jejbji 


ji...j6.'mi...m6 


The Feynman rules are the same than in the previous case, and the state sum amplitude associated to a graph is 
given by 


(-A)" ^ 1 = (-A)" ^ l[{2j, + 1). (31) 

fi.T { 3 } ® 

As before, this is interpreted as the amplitude associated to a 3D triangulation. This is a variant of the Ponzano-Regge 
model, with 1 as tetrahedron amplitude, instead of the 6 j—symbol. 

This partition function obviously diverges and has to be regularized. As before, an elegant way to regularize this 
model is to consider the original model to be defined on the non-commutative sphere instead of the 2-sphere. The 
algebras of function on Sff is just the space of {N -I-1) x {N + 1) matrices. This space carries naturally a representation 
of SU(2); if we decompose it in terms of spherical harmonics it can be written as a sum over spin j representation 
where j < N (see appendix B). The amplitude associated to a graph becomes 

•IPI - ^ E^^ n(2> + 0 = ^[(« + 1)7* (32) 

= ^ ( y)(# tetrahedra) (# edges)^ 2333 

sym[F] ^ 

which is finite (we defined fi = (A^-|-l)^). This amplitude is the amplitude (1) obtained in 3D dynamical triangulations. 
However, we have to notice that to correctly interpret this GFT as defining a sum over triangulations of the 3D 
dynamical triangulation model, we have to set A negative, namely —A ~ exp(—A — ^), see (3). 


C. Other group field theories 

The two previous group field theories were the first ones introduced in the literature [5, 8 ]. There has been since 
then a resurgence of these models when it appeared that the Barrett-Crane amplitude, which was a candidate for a 
4D state sum model of quantum gravity, could be interpreted as a Feynman graph amplitude of a group field theory 
[ 11 , 18]. 
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Since then, different modifications of the original 4D model have been proposed first by Perez and Rovelli [14] (note 
also the generalization in [19]). The purpose of these modified models was to improve the convergence properties of 
the Feynman graph amplitudes. We give here the 3D analogs of these models. 

The general group field theory we will consider is defined on SO(3)^ for real symmetric fields, with the action 

<S'[4>] = \ j dgidg2dgz [Pkin^{gi,g2,g3)f 

+ dgidg2dg3dg4dg5dge [Ppot'^igi, g 2 , gsMPpot'^igs, g^, gi)][Ppot‘^ig4, g 2 , g6)][Ppot‘^ig6, g5, gi)], (34) 

where Pkin, Ppot are operators acting on the space of fields and implementing some invariance property. For instance 
the Boulatov model is recovered for Pkin = Ppot = Psois)- Where Psoi3)^{gi, g 2 , ga) = !so( 3 )^^ 9 ig,g 2 g,g 3 g)dg 
is the projector on the space of S'0(3) invariant fields. We denote by Pso( 2 p projector on the sphere, i-e 
-^io(2)3^(ffi>52,53) = Jso( 2 ) 3 ^i 9 ihi,g 2 h 2 ,g 3 h 3 )dhidh 2 dh 3 . The Ambjorn model is recovered for Pkin = Ppot = 
PgO{ 2 p- consider new group field theories that gives 3 dimensional analogs of the models introduced by Perez 

and Rovelli, if we use 


Pkin — PsO{3)^ 

Ppot = {PsO{3)PsO{ 2)^Y Pkin- 

As the kinetic term is the same than in the Boulatov model, the Fourier coefficient will be the same, 
applying the projectors in Ppot leads to the following state sum amplitude 


(35) 

(36) 

However, 


“^F]=En(2j/ + l)ni f, % 2 ^e^’(ji,j2,j3)0"(j3,j5,j4)0"(j2,j4,j6)0"(j-l,j5,j6) 
{ 3 } f ^ ^ 


(37) 


where 0 is given by 


0(ji, 72 , 43 ) = \Civ^ = dg D^()(g)Dj'(,(g)Dj'"o(g). 


(38) 


One can use the following majoration of 0 


0(71,72,73) = 


\ ^ r^ji32j3 

— / . ^0m2m3'~^0m2m3 


(39) 


This last term in the RHS is just l/dji- Applying this uniformly for the three spins one gets the bound 


0(71,72,73) < 


1 

{djidj^dj^y/3' 


(40) 


III. NON-PERTURBATIVE RESUMMATION OF AN ASYMPTOTIC SERIES 

In this section we review well known facts on divergent series in quantum theory. We recall the physical origin of 
these divergences, and how to deal with them from the perturbative and non-perturbative point of view. In particular, 
we explain how to obtain a non-perturbative resummation of a divergent perturbative series. We illustrate these ideas 
on a one dimensional (fY model and show the crucial role played by the instantons in this framework (see [13] for a 
general overview of the large order behavior in quantum field theory.) 

It is known since Dyson [6] that physical quantities in interacting quantum theory are non analytic functions of the 
coupling constant. This is easily related to the fact that for negative values of the coupling constant, the potential fails 
to be bounded from below and the theory becomes unstable. The consequence is that the perturbative expansion of 
a physical quantity in terms of the coupling constant is a divergent series. A simple example is given by the following 
action 


Sxic^) = 



(41) 
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FIG. 3: Contour for the analytic continuation of Z{\) 


which is unbounded from below for negative A. The partition function 


/ + 00 
-OO 


can be perturbatively evaluated by first expanding the integrand in terms of A, then inverting the integration and the 
summation and integrating each term 


/ +00 1 

d<j) = 22"+5r(2n + -). 

-OO 2 


The exchange of integration and summation is in fact invalid, as a result the perturbative expansion of Z{X) 

\n 1 

Z(A)=^^/2(-l^—r(2n+-), (44) 

n 

is a divergent series with the n-th term growing like C(—l)"4"n!. 

In general, the asymptotic behavior of perturbative divergent series can be analyzed by considering classical non- 
perturbative objects of the theory called instantons. The instantons are complex valued solutions of the equations 
of motion with finite action. For instance considering the action (41), its equations of motion have two non-trivial 
complex solutions given hj (j) = ii/VX. The value of the action for these solutions is 5(0) = — These quantities 
are expressed using the inverse of the coupling constant, showing the non-perturbative nature of these objects. Let 
us see how they play a role in the large order behavior of perturbative expansions. 

Consider the previous partition function (42). First, we need to show that it can be analytically continued in the 
complex plane, with a cut along the real negative axis. Let A = be a complex number. We define Z(\) by 

Z{X) = f d(t) (45) 

First Z{\) is well defined for all \9\ < tt, second it is an analytic continuation of Z{X). This is clear since for \9\ < 
the integral along is equal to the integral along R plus the circular contours which link these lines at infinity 

(see figure 3). Writing z = i?e*“ on these circular contours, the integral along them is equal to 


lim 

R —»-±oo 




which goes to zero when R ±oo, provided that |0| < 7r/2. Thus the contour of integration can be shifted from 
g-ie/4]^ to M and the proposed expression defines the analytic continuation of Z{X), which can be extended to the 
whole complex plane, with a cut on the real negative axis. We can then define the discontinuity along this cut axis by 

D{-uj) = Z{e^^Lo) - Z(e-^^c^) (47) 
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FIG. 4: a) Original contour C for defining b) Deformed contour C'to apply the saddle point method. 


for a; > 0. Thus D is given by the integral 


D{-uj) = [ d<j) (48) 

Jc 

where C is the contour pictured in figure 4. Now Z is related to this discontinuity by the dispersion 

relation 


Z{\) 


ziTT Jq ct; + a 


We can expand this dispersion relation 


Z(A) = ^a„A" 

n 


with 


an 


(- 1 )" 

2f7r 



^n+l 


(49) 


(50) 


(51) 


We want to understand the large n behavior of the coefficients a„ of the series. For large n, the integrand in the 
definition of these coefficients is sharply peaked around z = 0. This shows that the large n behavior of a„ depends 
on the value of D{—uj) for small lo. We can change the variable of integration ip = -/oj<p. 


D{—u;) = [ d(p e 

Jc 



(52) 


The value of the integral for small to can be obtained by saddle point method. The saddle points equation is 
—ip + ip^ = 0, whose non-trivial solutions are ip = ±1. We see that these solutions corresponds to instantons solutions, 
since if we rescale ip> by iy/cJ, the saddle point equation becomes the equation of motion of the theory and the solutions 
becomes the instantons. The contour C can be deformed into C (see figure 4) which contains these saddle points. 
Expanding around these points and integrating, one gets 

D{-uj) oc (53) 


which is the evaluation of the action for the instantons solutions. One gets 


C+OO 

C(-l)” / diV 

Jo 


1 


U! 


,n+l 


(-1)” 7 


(54) 


which is in agreement with the asymptotic given for the exact computation (44). The value (1/4) " arises from the 
evaluation — l/4z of the action on instantons solutions. This computation shows clearly how the instantons govern 
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the behavior of the asymptotic series. We will now see how they are related to the way of extracting perturbative and 
non-perturbative information from the divergent expansions. 

Even if interacting quantum theories lead to perturbative expansions which are divergent series one can use these 
series as approximation of the non-perturbative quantity by cutting them at an appropriate order. Consider a divergent 
series 

N-l 

f{z)='^anz"‘ + RN{z) (55) 

n—0 

with a rest bounded by 

|i?iv(^)| < (56) 

(as seen in the previous example, S arise from the integration (54) and is related to the evaluation of the action for first 
instantons solutions.) One can use the series truncated at order N as an approximation of the exact non-perturbative 
quantity f{z), the difference between them being given by the rest Rn { z ). One can study the variations of Rn { z ) 
with N (at fixed z) and show that the rest first decreases as N increases, reaches its minimum around N ~ S/z, then 
increases and diverge to infinity. This means that the best estimation of f{z) by its truncated asymptotic series is 
obtained at order N ^ S/z. One can define the accuracy of the expansion by 

e{z) = Min^v Rn { z ). (57) 

Evaluating the bound (56) for N ^ S/z, one sees that e{z) ~ e~^l^. 

One thus know how to extract the best accurate physical information from a divergent perturbative series. However, 
one would like to be able to extract exact finite non-perturbative information from the divergent perturbative expan¬ 
sions. In order to get non-perturbative quantities from perturbative expansions, one has to reconstruct the original 
function Z{X) from an asymptotic series. However, it is easy to see that the issue of unicity of such reconstruction is 
not a trivial problem. For instance, if f{z) as an asymptotic series ^ akZ^, then f{z) -I- has the same one. This 

is not a surprise since one expects an ambiguity arising from the fact that the accuracy function is not zero, which 
means that the perturbative expansion does not contain all the information necessary to reconstruct the original 
function. Thus without additional conditions on /, it remains this ambiguity, there is no way to reconstruct uniquely 
the function from the asymptotic series. If one can find conditions on / to uniquely reconstruct it from its asymptotic 
series, the asymptotic series is said to be uniquely Borel summable. 

In that case, the unique function is recovered for instance by the Borel transform procedure. Let us consider an 
asymptotic series Assuming that the series 

converges in some circle |t| < r, and that its analytic continuation admits no singularity along the positive axis, one 
can consider the function 

^- 1-00 

f{z) = / e~*B{zt)dt, (59) 

Jo 

whose perturbative expansion gives the asymptotic series 'Y/akZ^. The possible ambiguity in the reconstruction is 
manifest in the Borel transform procedure in the fact that the integral (59) can not be unambiguously performed if 
the Borel transform B{t) admits poles along the positive axis. In this case we will have to specify how we go around 
the poles and different prescription will differ by terms of the type . 

The Sokal-Nevalinna’s improvement of Watson theorem [20] gives specific criteria on / in order for its asymptotic 
series to be uniquely Borel summable. 

Theorem 1 (Sokal-Nevalinna) Let f be analytic in the circle Cr = {z|3?(z“^) > R~^} (see figure 5) and having 
an asymptotic expansion 

N-l 

f{z)='^akz'" + RN{z) (60) 

with 

|i?7v(2)| <AS-^r{N + b)\zf 

then the asymptotic series is uniquely Borel summable and f{z) can be uniquely reconstructed from it. 


( 61 ) 
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FIG. 5: Domain Cr of analyticity for /, and the corresponding domain of analyticity for its Borel transform Bf, which contains 
an open neighbourhood of the real positive axis. 


The proof of the theorem can be outlined as follows. Suppose we have two functions /i and /2 having the desired 
asymptotic expansion and satisfying the hypothesis of the theorem. The difference between these two functions is 
given by their rests and one can bound these difference using (61). Our aim is to prove that this difference is null. 
Using the minimum on TV, it can be bounded by the accuracy of the asymptotic expansion 

e{z) ~ (62) 

Any such function analytic in the region Cr is in fact null. For instance the function is not analytic in Cr : 

let us consider the path z(t) = (which is in Cr for t < R). This path approaches z = 0 as t = 0 where as 

91(l/z(t)) = 2sint/t ^ 2 so I exp“'^/^*^*) | doesn’t approach 0. From the point of view of the reconstruction via the 
Borel transform, the analyticity in Cr ensures the analyticity of the Borel transform on a strip like region around the 
real positive axis (see figure 5). Thus there is no ambiguity in the choice of the contour along the real positive axis. 

Coming back to our example, with these criteria one can prove that the partition function (42) has an asymptotic 
expansion which is Borel summable. In fact we have already prove the criteria in our study of the behavior of the 
perturbative expansion : we have an analytic continuation in the cut complex plane which proves the first criterion 
and an asymptotic expression of the coefficients of the perturbative expansion which gives the desired bound for the 
second criterion. However, we will give now independent proofs of these criteria which are simpler and can be adapted 
to the GFT case. 

Concerning the analyticity in the domain Cr, if consider A G C, in the region 91(A) > 0, we have 


Thus we get the bound 




= e 


< 1 . 



(63) 


(64) 


which ensures the convergence of the integral for 9t(A) > 0. In particular, Z{\) is analytic in a region like the Cr 
needed in the first criterion. Concerning the second criterion, we start from this usual expression for the Taylor rest 


Rn{X) 


A^ 

(iV-1)! 




We use the same bound (64) for A G Cr to obtain 


d^Z 

dX^ 



< 

< 





4^ 


V2 r(2iv + i) 


This gives for the Taylor rest the bound 

\N 1 

\RNiX)\ < j^V2Ti2N +-), 


(65) 

( 66 ) 

(67) 

( 68 ) 

( 69 ) 
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Due to the multiplicative relation on the F-function T{2x) = ^/^r(a;)r(a; + 1/2), one can find a bound 

^ C'^^r(iv + b), (70) 

which leads a bound of the required form for the rest. The partition function satisfies the two Sokal criteria, hence 
its perturbative expansion is uniquely Borel summable. 


IV. MODIFICATION OF THE GROUP FIELD THEORY MODEL 


In this part we present a new type of group field models, related to the previous ones, and obtained by changing 
the potential term in the action. To determine this potential term, we are guided by the following ideas. First of all, 
the previous study of a simple quantum mechanical case has shown the crucial role played by the positivity of the 
potential in the proofs of Borel-summability. The tetrahedral potential of action (9) has to be modified since it is not 
positive. Moreover, it is natural from renormalisation group arguments to include all the potential terms of identical 
degree, which are allowed by the symmetries of the model (see appendix C for a comment on the divergent graphs 
which a priori need to be renormalized). From the point of view of the Feynman graph expansion, this is mirrored in 
the fact that the tetrahedron (which appears in the interpretation of the Feynman graphs as triangulations) is not the 
only natural three dimensional object constructible from 4 triangles. One can construct another object, the ’’pillow” 
(see figure 6). 



FIG. 6: The tetrahedron and the pillow : the two objects constructible from 4 triangles 

Our strategy is thus to include the corresponding interaction term in the GFT. We consider a general interaction 
term 


+ S 


which corresponds to a total action 

= ^ [ $(31,52,53)^ 

2 Jg^ 

+ 77 / ^{91,92,93)^(93,95,94)^(94,95, 96 )^( 96 , 92 ,91) 

4! Jgo L 

+ <5 $(51,52,33)^(33,55,54)^(34,32, 36 )^( 56 , 35 ,51) ( 71 ) 

Recall now that we are looking for a positive potential. We will determine the values of 6 giving a positive potential. 
The strategy is to rewrite the total potential as a scalar product involving a linear operator whose eigenvalues are 
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FIG. 7: A pillow is obtained by the gluing of two tetrahedra along two common faces, introducing an additional internal edge 
k 


under control. Let us define the object 


The pillow term is found to be 


A{9ii92,9b,9i) 


d 9 ^(51,52,5)^(5,55,54) 


(72) 



- 4 ( 51 ,52, 55 ,54)^4(51,52, 55 ,54) 


{A, A) 


while the tetrahedron term can be written as 


(73) 



Mai, 92, 55, 9 i)[T - 4 ](54,52, 55 ,51) 


{AMM 


(74) 


where T is the just the linear operator acting on A 

T :- 4 (^ 1 , 32 , 55 , 54 ) ^-4(34,32,55,51)- ( 75 ) 

This operator T is a permutation, thus = Id and its eigenvalues are simply ±1. The total interaction term we 
defined is thus 


{A, [Id + 5T]A) (76) 

whose eigenvalues are 1 ± (5. We obtained a positive potential for |(5| < 1. This complete the construction of the 
modified potential term. 

At first sight, one could say that such a modification of the potential term will lead to a drastic modification 
of the models, first because we will sum over different kind of triangulations, second because the weights will be 
different. We will see that the modification is not so drastic. We can first observe that geometrically, a pillow is no 
more than the gluing of two tetrahedra, along two common faces (see figure 7). This allows us to relate in unique 
manner each ’generalized’ triangulation (made of pillows and tetrahedra) to a ’real’ triangulation (made only of 
tetrahedra) by replacing each pillow by a pair of tetrahedra glued along two common faces. The total sum generated 
by the perturbative expansion can be reorganized by putting together the weights of all the generalized triangulations 
related to a given real triangulation, and summing only over the real triangulations. The real triangulations containing 
pairs of tetrahedra glued along two common faces are called ’irregular’, while the other ones are called ’regular’. We 
can see that this reorganization concerns only the irregular triangulations. 


Our results are the following. We consider the modifications of the Ambjorn and Boulatov GFT given by the action 
(71). In both cases the perturbative expansion of the partition function is uniquely Borel summable. In the case of 
the modified Ambjorn model, the perturbative expansion is written 


z{\M= E 

AeTres 


Sym[A] 




k 


( 77 ) 
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where /i = (A^ +1)^, N being the cut-ofF. The first sum is over regular triangulations, while the second is over irregular 
ones. (resp. Ni) is the number of tetrahedra (resp. edges) in A, k denotes the number of pairs of tetrahedra glued 
along two common faces in the irregular triangulations. In the modified Boulatov model, the perturbative expansion 
can be written 


where n is the number of tetrahedra in the triangulation A, k is the number of pairs of tetrahedra glued along 
two common faces in an irregular triangulation and Zp}i[A] is the regularized Ponzano-Regge amplitude, for the 
triangulation A. These two formulas have been written for the parameter 6 = —1 in the potential (see 71), but we 
will give a general expression valid for |(I| < 1 in the proof in the next section. 

Several remarks can be made. First, despite the modification made on the GFT models, we still obtain a perturbative 
expansion which can be written as a sum over triangulations. Moreover, for all regular triangulations which do not 
contain any pair of tetrahedra glued along two common faces, these weights are the same than the original models. 
The weight are modified for the irregular triangulations and can even be cancelled by an appropriate choice of the 
coupling constant, namely A = 1/2 in the Boulatov model or /xA = 1/2 in the dynamical triangulation model. In 
addition, our modification lead to another surprising consequence. The modification of the potential we choose leads 
to the fact that the factor depending on the coupling constant is now A" instead of (—A)", as it was the case for 
the originals models. Thus the relevant regime is for positive values of A, namely the physical region of convergence, 
where the GFT is well-defined. Our choice of positive potential allows us not only to prove the Borel summability of 
the model, but also to have a well-defined GFT for the values of the coupling constant which are physically relevant. 


V. PROOF OF THE BOREL-SUMMABILITY OF THE PERTURBATIVE EXPANSION 

In this part we prove the main result, namely the unique Borel summability of the two modified models, and the 
fact that the resulting perturbative expansions are rewritten as slight modifications of the original models. The proofs 
of Borel summability are closed to those given in the simple quantum mechanical case at the end of section III. 


A. Sokal criteria for the modified models 


The analyticity of the partition function for 1H(A) > 0 follows directly from the positivity of the potential, using 
the same bound than in the case, eq.(64). This proves the first criterion. 

We consider now the Fourier expansion of these models. For the 3D matrix model, we find 




m|2 

2 ^ '--3 I ' 4! 

J,fn 


A— mi—m2 — jms —77157714 J—771477157716 J — 77 l 677 l 2771 i 
^jlj2j3 ^j3j5j4 ^363231 


_1_T 7711—7712 —7713 47713 —77I577I4 4—771477127716 4 —77l677l5771i 

^jlj2j3 ^j33534 ^j4j2j6 ^jejsjl 


For the Boulatov model, we find 


^=pE 1^7 


4! 


E 


(it+771i) 


3,n 


^1,...J6,Till,...7716 


—7711 —771 2 —7713 /I 7713 —7715 7714 ^ —77I477I577I6 /I —77l677l2771i ^ O O O ^ 33^36 

^313233 ^333534 ^34h3G ^363231 ^3 i323q^34^33^343^3q . 


’ 2j3 + 1 


_l_r 4—7711 —7712 —7713 477I3 — 77I577I4 4—771477127716 4—771677157711 J 1 J3 

^313233 ^333534 ^34323G | jg 


(79) 


(80) 


where Oabc = 1 if a, 5, c satisfies the triangular inequalities, 0 if not. 

In order to clarify the notations and give a uniform proof for the two cases, we separate the coefficients into real and 
imaginary parts. All the real coefficients we obtained are not independent due to the relations ( 16 , 17 ). We denote by 
{A/, I = 1 ...C?} a set of independent coefficients {d is related to the cutoff N and behaves roughly as V®, as it stands 
for the six indices (ji,J 2 ,J 3 , wi,m 2 ,TO 3 ), up to symmetry relations). The action expressed in Fourier modes can be 
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rewritten as 

+ AiAjM^^^^AkAl (81) 

I IJKL 

with a four indices matrix representing the potential, whose positivity is translated into 

AiAjM^-^^^AkAl > 0 WAi (82) 

IJKL 

The partition function is 

Z{X)= f exp A ^ 

L / J I / IJKL } 

We want to consider the asymptotic expansion of Z{X) 

n—1 

z{X) = Y^ pXP + Rn{X) (84) 

p—0 

We use the expression (65) for the Taylor rest. Consider the expression for obtained by deriving n times the 
original expression (83) 

(At) = f []^dA/] exp ~ ^ E AiAiM^-’^^AkAl Y AiAjM^-^^^ArA^ (85) 

One look for a bound of the module of the integrand. The real part of the potential is positive for A € Cr, thus 
exp E AiAjM^-^^^AkAl < exp -^E^? 

Thus one obtains the following bound 

r ~\ 

/jn 7 r ___ -I __ -| 

-^(At) < / [J^dA/Jexp it E AiAjM^^^AkAl (87) 

I L / J IJKL 

Let us denote 

S = ( 88 ) 

One gets the bound 

jn 7 on r i / \ 

dvr(V) <m 

The RHS integral is symmetric over the Aj, thus one can write it as d! times the integral over the domain Ai > A 2 > 

... > Ad (recall that d denotes the total number of Fourier modes and depends polynomially on the cut-off N). In 

that case, 

Y. \AiAjAkAl\ < (90) 

IJKL 

Then we can decouple the integrals over the Aj on LHS of (89) and extend the integral for J ^ 1 in order to get d — 1 
gaussian integrals. The integral over Ai gives a F—function 






( 91 ) 
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and 


d^Z 

dA" 


{\t) 


< 


d-1 
2 d! TT 2 



r(2n+i) 


(92) 


Coming back to the expression of the Taylor rest (65), one gets a bound 


|i?„(A)| <C 




(93) 


which is of the required form, with constants depending on d, hence on the cut-off N. The perturbative expansion of 
our model satisfies the two Sokal criteria and thus is uniquely Borel summable. 


B. Reorganization of the perturbative expansion 

We now analyze the consequences of the introduction of the pillow interaction term on the expression of the sum 
over triangulations. The perturbative expansions of modified GFT models give us a sum over triangulations made 
with pillows and tetrahedra. However, we have seen that geometrically a pillow is obtained by gluing together two 
tetrahedra along two common faces. 

In the matrix model case, the amplitude associated to a pillow is 

Apiiiow = —A (94) 

while the amplitude for two tetrahedra glued along two common faces (see figure 7) is 

A^tetra = ^ ( 2 fc + 1 ) = {5\f [N + if (95) 

k<N 

Thus, if we compare the total amplitude for a (generalized) triangulation A containing k pillows to the related real 
triangulation A obtained by replacing all the pillows by pairs of tetrahedra, we obtain (recall fi = {N + 1)^) 

The one-half factor comes from the fact that the symmetry factor for two tetrahedra is twice the symmetry factor 
for one pillow. Consider a real triangulation A made only with A 3 tetrahedra and Ni edges, and such that 2k of its 
tetrahedra are glued by pair along two of their common faces. There are several generalized triangulations related to 
it : (binomial coefficients) triangulations with one pair of such tetrahedra replaced by a pillow, triangulations 

with two such pairs replaced by pillow etc...Counting together all their weights we obtain 

Z[A] = {-^X)^^(l + {-26‘^X^i)-^Cl + {-2^‘^Xfi)-‘^Cl + ... + {-2^‘^X^l)-^C^'^fi^^ = {-SX)^fl-i26‘^Xf^)-f'^^l^^ (97) 

This shows that one can reorganize our sum as a sum over real triangulations of the amplitude {—6X)^^ii^^, except 
the fact that an irregular triangulation containing k pairs of tetrahedra glued together along two common faces has 
a weight modified by a factor (1 — {2S‘^Xfi)~^f. Which gives 

2(A..) = E + s)* ^ 

Recall that the positivity of the potential is valid for |(5| < 1. For instance for <5 = —1, we recover the result (77) 
announced in section IV. 

The same reorganization can be performed on the perturbative expansion of the modified Boulatov model. The 
amplitude associated to a pillow is 


A 


pillow 


-X 


J 6 

2J3 + 1 


(99) 
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The amplitude associated to the gluing of two tetrahedra along two common faces (see figure 7) is 


A2tetra — + 1 ) 


Ji J2 k 
ji 35 J3 


ji 32 k 
ji 35 je 


= (^A)^ 


2J3 + 1 


( 100 ) 


due to the orthogonality relation. As before, this means that a triangulation containing a pillow has an amplitude 
which is (1 — 2d^X)~^ times the amplitude of the same triangulation, with the pillow replaced by two tetrahedra glued 
along two common faces. We reorganize the sum in the same way than before and obtain 


m = V 


The result (78) given in section IV is recovered for 5 = — 1. 


VI. CONCLUSION 

In this paper we have studied the possibility of defining the sum over topologies of three dimensional gravity in 
the context of dynamical triangulation or spin foam models. We have provided an example of a group field theory 
model which is uniquely Borel summable in the physical regime of the theory. In order to do so we have used the 
fact that three dimensional triangulations are generated by a group field theory with a quartic potential and we 
have shown that this potential can be naturally modified in order to be positive. Therefore, the message we wanted 
to convey here is that it is possible to resum the sum over triangulation if one accepts to take for serious the non 
perturbative information provided to us by the GFT, namely the instantonic solutions. One can still question what 
is the physical content of this non perturbative information. This is a direction of investigation we are following. 
Finally, 4 dimensional triangulations are generated by quintic potentials and we may think that there is no natural 
modification of this potential leading to a uniquely resummable theory. This is an issue that deserves to be studied 
and could mean that we need more non-perturbative input in this case in order to give a meaning to the resummed 
series. 
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APPENDIX A: REPRESENTATIONS OF SO(3) AND RACAH-WIGNER SYMBOLS 

The irreducible unitary inequivalent representations of SU(2) are labelled by a half-integer j. The representations of 
SO(3) are just the integer spin representation of SU(2). The space of representation is denoted Vj and has dimension 
dimj = 2j -I- I. All the representations are self-dual V* = Vj. The matrix of representations are denoted D^^{g). We 
have the relation 

DLuig) = (Al) 

Given three representations it exists only one (up to normalization) intertwiner i : ^ Vj^ G Vj^ C. 

One can express it by using the Wigner 3j-symbol whose normalization is 

E44 = 1 (A2) 


The Racah-Wigner 6j symbol is defined by 


f ji 32 js 
1 ji J5 J6 


' (_-\'\3i+jii+j6+'mi+m^+mef~<jij2j3 

■*-/ — —17147712 ^m4 — m5m3 


It satisfies the orthogonality relation 


E(2J3 + 1)(2J6 + 1) 

On 


jl J2 J3 
ji J5 J6 


jl J2 j's 

ji J5 J6 


= (i 


Jsj's ^010203 ^J3J4J5 


(A3) 


(A4) 
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where 0abc = 1 if a, 6, c satisfy the triangular inequalities, 0 in the contrary case. The fact that the integral of a 
matrix representation gives a projection on the space of invariant vectors leads to the following relations (using the 
normalized Haar measure) 


and 



dgDi,^{g)D^^,n,{g) 


_ ^33 ^ , 

dj 


(A5) 




02 

m2n2 


(g)Di^ „ (o) = ClCi. 

Vc'/ m^nz\^J m n 


(A6) 


APPENDIX B: THE NON-COMMUTATIVE 2-SPHERE 


A non-commutative space A is a space which is defined by specifying the (non-commutative) algebra of square 
integrable functions on it In this part we recall some facts about the construction of the non-commutative 

sphere. The sphere S'^ can be described by the algebra A basis of this algebra is given by polynomes of 

such that 

x\+ x\+ x\ = \ (Bl) 

Another possible basis is given the spherical harmonics F*™ and arises from the decomposition 

+ 00 

L\S'^)=^Vi. (B2) 

1=1 

where Vj is the spin j representation of SO(3). Let us construct the non-commutative sphere (for a deformation 
parameter 1/iV) by giving its algebra of functions, obtained by starting from the first basis described for the sphere. 
Let us consider Vjv, the spin N/2 representation of the algebra su(2); are the matrices of representation. Let 

us define 

(B3) 

where the Ji are the generators of su(2). Xi is in V_n. Due to the Lie algebra relations of su(2), one gets 

= 1 (B4) 

i 

[Xi,Xfi = (B5) 

The first relation is the analogue of the relation (eq. Bl), and the second one says that the Xi are commuting in the 
large N limit. This facts lead us to see the Xi as non-commuting spherical coordinates, and to define the algebra of 
functions on the non-commutative sphere as the algebra of endomorphism of Vjv 

= EndiVn)'^V n(^Vn (B6) 

The right-hand side can be written as 


N 

(B7) 

1=0 

which shows that the parameter N is nothing more than the cut-off in the decomposition of the functions over 
the sphere into spherical harmonics. As the endomorphisms of VA are the non-commutative counterpart of the 
polynomials over the sphere, a natural question is to ask what are the non-commutative counterpart of the spherical 
harmonics basis. Indeed, each spherical harmonic for I < N can be represented as an endomorphism of Fw by 

l,m _ 2 2 

^ Hj ~ i j • 


F' 


(B8) 



20 


APPENDIX C: NOTE ON THE CONVERGENCE OF GRAPHS 


If we consider the amplitude associated to a graph F in the group field theory obtained by using p times the 
projection (see subsection IIC) , we get 

^[ri=i:ii‘i«ne"oi,A,A)n{(; <ci) 

{j} f e V ■> 

We have proven for the 0 the following bound 




1 


We consider also the bound on the 6j-symbol 


f Ji J2 js 

\ J4 jb je 


~ i-n& 




obtained from the orthogonality relation (A4). 
The amplitude (Cl) can be bounded 


4r]< 


li/l / 


(C2) 


(C3) 


(C4) 


where Ueif) and ny{f) are the number of edges and vertices lying in the face /. As these numbers are equal, we 
obtain the bound 


-4[r] < n 

/ 



(C5) 


One can see that there is no divergence associated to the sum over spins of face / if the number of vertices lying in 
this face is such that 


nv{f) > 


12 

2p+ 1' 


(C6) 


For instance in the case of p = 2, we obtain that the only graphs which are not proved to be convergent are for the 
faces such that n^(/) = 1, 2. The dual faces bounded by only one or two dual vertices correspond to edges bounded by 
only one or two tetrahedra in the triangulation. This observation justify that only these two types of graphs need to 
be renormalized. However, our bounds are not optimal and numerical simulations suggest that all the graphs actually 
converge for p = 2, while for p = 1, the graphs with nf(v) = 1,2 diverges. 
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